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Abstract 

Some exact analytical formulas are presented for the generalized Roche 
model of rotating star. The gravitational field of the central core is de- 
scribed by the model of two equal-mass point centers placed symmetri- 
cally at rotation axis with the pure imaginary z coordinates. The all 
basic parameters of the critical figure of the rotating massless envelope 
are presented in analytical form. The existence of the concave form of the 
uniformly rotating liquid is shown for a large enough angular velocity of 
the rotation. 



1 Introduction 

The classic Roche model describes the equilibrium figure of the rotating massless 
liquid (gas) envelope in the gravitational field of the point mass [1,2]. It is vifidely 
used in investigation of the structure of the rotating star with a large enough 
value of the effective polytropic index, when there is a strong concentration of a 
matter to the star's center and one may ignore the self-gravitation of the outer 
layers of the star as compared with the gravitation of the star's central core. In 
the opposite limit of rotation of the incompressible liquid with constant density, 
the Maclarain ellipsoidal models are used [1,2]. 

The next natural generalization of the classic Roche model would be the taking 
into account the effects of the non-spherical gravitational field of the star's core, 
still neglecting the self-gravitation of the (uniformly) rotating envelope. In this 
note we consider one of the possible approaches to the problem [3] . 



2 Basic Equations 

2.1 Gravitational potential of two fixed centers 

We use a common set of rectangular coordinates: x- and y-axes are in the equa- 
torial plane, z-axis coincides with the axis of rotation, the center of coordinate 
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system (0, 0, 0) coincides with the star's center. 

Let two fixed points with equal masses mi = m2 — M/2 be placed at z- 
axis such that coordinates of the ^centers^ are (0,0, 2:1) and (0,0,Z2), where 
zi = ic, Z2 = —i c, = —1. 

Then gravitational potential has the axis of symmetry, z-axis, (axis of rotation), 
and the plane of symmetry, z = plane, (equatorial plane). For our problem 
it is sufficient to consider the gravitational potential in (a;, z)-plane (meridional 
plane) 

V{x,z)^ ^ {- + -), 1 

rl= x^ + {z-icf, rl=x^ + [z + i c)^, (2) 

where G is the gravitational constant, M is the star mass, ri and r2 are distances 
from the point (a;, z) to fixed centers, c is the (positive) constant of the two fixed 
centers problem describing the deviation of the gravitational potential from the 
spherical symmetry. Hereafter all formulas reduce to the classical Roche model 
case at the limit c ^ 0. 



2.2 Reality of gravitational potential 

Gravitational potential (|1I2|I is of course real (while ri and r2 not), and in 
order to show it we introduce, instead of rectangular coordinates {x, z), two real 
positive variables A and /i defined by the following relations: 

ri=c{\-i^i), r2 = c(A + i/i), (3) 
ri + r2 r2 - ri 



2 c ' 2 c 



(4) 



X 



From Equations l^-© we get 

.r,. N GM A 

^('''^)~A^' 

which is real, QED. 

Note that (A,/x) coordinates are sometimes called (oblate or prolate) confocal 
ellipsoidal coordinates, or simply ellipsoidal coordinates, or even elliptic coor- 
dinates, see e.g. [4,5]. Also notes that authors use different notations for these 
coordinates. At (a;, z)— plane, the curves A=const are confocal ellipses with fo- 
cuses at points (a; = 0, z = +c) and (a; = 0, z = — c) while the curves /i=const 
being normal to curves A=const are hyperbolas. 
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2.3 Total potential 

Potential of rotational forces is 

U{x) = ^u;'x^ (7) 

where oj =const is the angular velocity of rotation. The total potential is a sum 
of gravitational and rotational potentials 

<i>{x,z) = Vix,z) + U{x). (8) 

Again we note that the total potential has an axis and a plane of symmetry 
though we consider only {x, 2;)— plane (meridian plane) and more specifically, 
first quadrant of {x, 2:)— plane with positive x and z. 



3 Arbitrary rotation in the classic Roche prob- 
lem 

The form of the equilibrium rotating masslcss envelope in the CRM depends on 
degree of rotation: non-rotating envelope has spherical form, for slow rotation 
the figure is spheroid (figure of rotation with meridional section as ellipse), for 
larger rotation the figure is more and more oblate (and not spheroid), imtil 
the critical figure is attained, when at the equator of the figure, the centrifugal 
force is equal to gravitational force. Note that critical figure has a cusp at 
equator. For arbitrary value of angular velocity u the total potential at the 
equator (where r = a; = a) is 

GM 1 2 2 
$c = — + « • (9) 

The equilibrium figure is defined by the equation 

GM 1 2 2 GM 1 2 2 nnA 

h -io^x^ = h -OJ a^, (10) 

r 2 a 2 

from here we find the polar radius {x = 0,r = z = b) 

GM GM 1 



6 a ^r^' (11) 

for the ellipticity / = 1 — a/6 of the equilibrium figure and centrifugal- to- 
gravitational force ratio at equator m = ui^a/{GM/a'^) we have the general 
expression 

f 1 

- = 7^ • (12) 

m 2 + m 

We see that, in the CRM, the ratio //m varies from maximal value of 1/2 (for 
slow rotation, m <C 1), to minimal value of 1/3 (for maximal rotation, m = 1). 
Note that sometimes authors use in m centrifugal force at equator and mean 
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gravitational force which is OK in the first approximation in cj^ but has no 
much sense at strong rotation. 

For comparison, for Maclaurin spheroids, the //m ratio varies from 5/4 to 1 for 
m varying from to 1, see Fig. ^ 




Figure 1: Elhpticity of rotating fluid in Maclaurin and Roche models: abscissae 
is m, centrifugal-to-gravitational force ratio at equator of the figure, ordinate is 
//m, where / = 1 — 6/a is cllipticity, and a and b are equatorial and polar radii. 



4 Critical figure in the GRM: 
Polar and equatorial radii 

Let us write down the gravitational potential (1,2) and the total potential (jSJ 
on a;-axis, where z = 0, ri — r2 — \/ — c^, a; > c > 

C M 

V{x)^ <i>{x)^V{x) + U{x). (13) 

Vx^ — 

The critical surface is defined by a condition of equality of centrifugal and gravi- 
tational forces, -dV{x)/dx = dU{x) /dx, or by a minimum of the total potential, 
d^{x)/dx = 0, at some x — a. This gives the equation for the equatorial radius 

(14) 



(a2 - c^f^ 
4.1 Equatorial radius 

From Equation H14|l. the equatorial radius of the critical surface is 



+ c2. (15) 
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We note that the equatorial radius of the rotating star with the same mass M 
and angular velocity lo is larger than the equatorial radius oq in the classical 
Roche model 

ao = — ^ . (16) 



4.2 Polar radius 

From equations l(7jl- p5(l . we write the total potential at the critical surface as 

GM 3a' -2c' cj^ 
$0 = ^ =— = — {3a'-2c'). (17) 

Y (a' — c') 

On rotational, z, axis, where x — 0, ri = z — ci, r2 = z + ci, z > c > , we 
have only the gravitational potential (because the rotational potential vanishes 
at z-axis) 

V{z)^GM^^. (18) 

The polar radius of the critical figure is defined by the condition $o — ^iz), 
and from equations (|17|l and H18I) , we find the polar radius of the critical figure 



{a' ~ c2)3/2 



Note that we choose plus sign before radical in (|19|l such that b > c, also we 
note that usually the value of the constant c of the two fixed centers model is 
much less than dimensions of the figure (star or planet). 

From (|14|l we have unequality for the equatorial radius a > c, additionally from 
condition b > c, {A > c), and (|19|l . we get the lower boundary for the equatorial 
radius 



ai 

a > ai, — = 
c 



^\ 1/3 / ^\ 1/3 

73-V5\ /73 + V5 



Numerically, ai/c = 3.26092 [5]. 



4.3 Two limits 

At limit of small deviations from the classical Roche model, a ^ c, b ^ c, we 
have for equatorial, a, and polar, b, radii and their ratio b/a 

f I c'\ [gmV'^ 
31 c' 

^0 



6 = 6o(l--^ , 5o = -ao, (22) 
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37 
18 



1025 
216 



(23) 

We note that comparing with CRM, in GRM, the equatorial radius a is larger, 
while the polar radius h and the polar-to-equatorial radius ratio hj a is less. That 
is in GRM, the critical figure of rotating envelope is more flatten (oblate). 
In another limit, at the equatorial radius a close to we have 



a 



- = ^+S, 0<S<^: 
c c c 



= 1 



GafcS 



-,1/2 



^ = ^fl-^ 
a ai \ ai 



(a2_c2)(3a2-2c2) 
6alSc 



^ ~ (a2-c2)(3a2-2c2) 



11/2 



At (5 = 0, a = fli, 

b 
a 



\5{5 + V5)) 5^5(5 + ^5) 



(24) 

(25) 
(26) 

(27) 



5 Isopotentials in GRM 

In general form the equation for critical figure can be written down using vari- 
ables A and ^ 



GM A 



1 



A2 



+ -l.^c^(1-HA^)(1-m') = *o, 



(28) 



where $o and a are given in H17|) and (|15|l . 

From (|17|l we have quadratic equation for fi^ as function of A. Further, in vari- 
ables X, z using Q we get parametric equation of the complex from (with A as 
parameter). In result, we present the Fig. 2 the equilibrium critical figures for 
several values of the equatorial radius a/c. With decreasing a (this correspond- 
ing to increasing w or increasing c, increasing flatteness of the gravitational 
field) the deviation of the GRM critical figures from CRM increases and for 
small values of a the figures become even concave while CRM figures are always 
convex. 



6 Concave figure 

Let us find the value a = ai when the equilibrium figure becomes concave first. 
Expanding (1,2) in the series at small a; ^ 1 and using (7,8,17), and requiring 
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dz/dx = a,t X = we have the algebraic equation of the 7th order for y = 02/0 



144-1512 y+6408i/2-141842/3+i7475/_ii553y5+3397 2/6-171 y7 ^ q. (39) 

similarly for corresponding polar radius 62 we have again algebraic equations of 
the same order {t = hijc) 

81 + 584 i + 1278 + 1500 + 882 1'^ + 195 - 45 - 19 1' = 0. (30) 

Numerically, we get a2/c = 4.00358, 62/c = 2.0376, and 62/012 = 0.508945. 

At a > 02 the critical figures in GRM are convex, at a < a2 the critical figures 

(and also inner isopotentials near the external surface!) are concave. 
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Figure 2: Critical figures in GRM for different values oi a/c. At a/c > 4.0 
equilibrium figures are convex, at 3.26< a/c <4.0 the equilibrium figures are 
convex near the polar axis. Also shown are the equilibrium figures in CGM for 
the same values of a. 



7 Post-classical approximation 

The analysis is more easy in the limit x c, z c (more exactly, r = (x^ + 
^2^1/2 rpj^g gravitational potential up to terms proportional to c^/r^ is 



V{r,z) = 



GM 



1 + 



3 C2 ^2 



2 ^2 2 r2 



(31) 
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Isopotentials V{r, z) — Vc corresponding to are ellipses with larger semi- 
axis a along a;-axis and minor semi-axis b along z-axis 



lc\ . c\ GM 

'o 

At c = we have the equilibrium figure of the CRM 

G M 



« = ^^0(1 + Tj-), b = ro(l - -), ro = — . (32) 



(33) 



or 

ro(x) = g _ J^^^ , < .T < a, a > ro > b. (34) 

If <f>c >> (slow rotation), the equilibrium figure is ellipse 

G M ,^ fcj2a;2^ 

$c 2 <f>c 

with larger semi-axes along the a;-axis (equatorial radius) and minor semi-axis 
along the z-axis (polar radius). 

Taking into account terms up to /a^ critical figure is 

r{x) = ro(x)(l -f /3cVao), < a; < a, 6 < r < a, (36) 

/3 = -^ + ^t'-3i'-i', i = aoAo, (37) 

where ro(x), ao,a, and b are given in Eqs (21-23, 34). 

Using (36,37) we may calculate the volume of the critical figure 

pa 

W^AttI z{x)xdx, z = Vr^-x^. (38) 
Jo 

In (38) both the integrand z{x) and the upper limit a can be expand in the 
series up to terms (? ja^ 



W = A-K 



ao 



zo(a;) a; dx, -l-z(ao) • ap • (a — oq) -f / z\{x')xdx 
h) 

2 2 2 I 2/2 all 



(39) 



Here the first term in W is the volume of critical figure in CRM, second term 
is proportional to (? ja^ and may me omitted, and third term is proportional to 
c2/a2.^ 

Both integrals can be expressed in terms of elementary functions, and we have 

r°-o fi 1 _ 7.2 4 / \ 

/ ZQ(x)xdx^al \ ^(4-x2)i/2a;da; = V3---Hln 6-3%/3 , (40) 

Jo Jo 3-x^ 3 V y 
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ra-0 

I zi{x) X dx 





16-27^3 
210 



(41) 



The volume of the critical figure is 



W = Ana] 



,1 [0.180372- 0.146502(cVao)] 



(42) 



By introducing the mean density of the critical figure p = M/W, we have in the 
post-classical approximation for parameter widely used in the theory of rotating 
configurations 



8 Conclusion 

We obtained all parameters of rotating configuration in the Roche model taking 
into account the effects of (small) non-sphericity of the gravitational field of the 
central body. As a rule any modification of the classical models leads at best to 
ugly non-elegant analytical or even pure numerical results. 
And Roche models are of the best examples of those classical models idealized 
and elegant both in set up and results. 
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0.360744 -0.293003 cVao- 



(43) 
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